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Abstract: Ising spins put onto a Barabasi-Albert scale- free network show 
an effective phase transition from ferromagnetism to paramagnetism upon 
heating, with an effective critical temperature increasing as the logarithm of 
the system size. Starting with all spins up and upon equilibration pinning 
the few most-connected spins down nucleates the phase with most of the 
spins down. 
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Networks with more complicated connectivities than periodic lattices have 
been investigated in detail recently. For example, the Barabasi-Albert net- 
work [|TJ is grown such that the probability of a new site to be connected to 
one of the already existing sites is proportional to the number of the previous 
connections to this already existing site: The rich get richer. Similar net- 
works have been investigated || to check if destruction of a few computers 
will split the percolating cluster of the internet (i.e. the set of all computers 
in the world connected directly or indirectly with each other). Networks ex- 
ist also in social models where vertices are individuals or organizations, and 
links correspond to social relationships between them B. However, as far 
as we know, all studies of the scale-free Barabasi-Albert networks considered 
only the topology and no interactions between linked vertices. 

Here we investigate the ordering phenomenon in this Barabasi-Albert 
network, if Ising spins are put onto the sites (vertices) of the network. We 
assume ferromagnetic coupling between linked spins and positive temper- 
ature T of the system. Such a magnet would show paramagnetism if the 
whole network is still connected but only weakly such that thermal fluctu- 
ations destroy the spontaneous magnetization. (Bose-Einstein condensation 
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on similar networks was already studied before J3J, as was the Ising model 
on small- world networks ||, and other Ising models are in preparation [[J.) 
In the social example one could identify exp(— const /T) with the probability 
that members of the same social group are not convinced to share the same 
opinion. (Spin variables Si = ±1 then correspond to two possible opinions 
of the group members on the same subject.) 

Thus we create a Barabasi-Albert network of iV sites added to an ini- 
tial core of m fully connected sites. Each of the N new nodes is connected 
to m randomly selected previous nodes. (We allow more than one connec- 
tion to the same site for the later added nodes; mostly we take m = 5.) 
Then we freeze this network structure, put an Ising spin Si = ±1 onto every 
site, with all spins up initially. Then with the standard heat bath Monte 
Carlo algorithm spins search for thermal equilibrium at temperature T (all 
temperatures are given in units of coupling constant over Boltzmann con- 
stant). Fig. la shows the resulting magnetization (averaged over the last half 
of 500 sweeps through the network) versus temperature; it seems to decrease 
exponentially with increasing temperature, until due to finite-size effects it 
oscillates about zero. This effective Curie temperature T C (N) is fitted in 
Fig.lb onto 2.6 ln(iV) - 3 for 5 < iV < 5, 000, 000. 

Figure 2 shows that in the ferromagnetic region the spins with few neig- 
bours flip up and down while those with many neighbours point up most of 
the time. 

Analogous to the appearance and spreading of new opinions in society |7J 
we now try to flip the spontaneous magnetization for T < T c by forcing the 
most-connected spin permanently down; then we pin in the same way the 
second-most-connected spin, and so on, all in time intervals of 50 iterations. 
We see that in general removal of the few leading spins having hundreds of 
neighbours is sufficient to flip the magnetization. (Since we do not apply any 
magnetic field and have no fixed boundary conditions, we expect the flipping 
of the magnetization to be a nucleation event, which would happen even if we 
flip only one randomly selected spin provided we wait long enough.) Fig. 3 
shows the magnetization versus time averaged over 100 samples. Higher 
temperatures require fewer leading spins to be pinned. 

Now we present a simple mean field theory for some of these effects. 
Let us consider the BA network with the characteristic constant m and the 
corresponding Ising model with the ferromagnetic coupling constant J = 1. 
The probability that a node has a degree k is given by P(k) ~ Ak~ 3 when 
k 3> m and P(k) = when k < m. For large networks (N — > oo) the 
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normalization constant equals A ~ 2m 2 . In the mean field approximation 
(MFA) we can simplify interactions among each group of spins with a fixed 
value of k by the effective field kM where M is the mean magnetization (per 
one spin). It follows that 



(/3 = 1/T). This is a transcendental algebraic equation for M(T). It is easy 
to find analytically the critical temperature T c that corresponds to the case 
when lhs and rhs are tangent at M — . Differentiating both sides over 
M and putting M = we get 1 = A j3k~ 2 dk; thus T c = 2m. The result 
can be written as T c =< k >~ 2m = the mean degree of the node. This is 
the typical MFA result. Fig.l shows it to be correct in order of magnitude 
and in its increase with increasing m. But the MFA does not describe the 
logarithmic size effect on T C (N), due perhaps to exponentially rare and small 
regions of densely connected spins 

Knowing the mean magnetization M(T) one can calculate the correlation 
fi(k) of local magnetization to the spin degree k. In fact in MFA we can 
write fi(k) = tanh[/3/cM(/3)]. It follows that for all temperatures T < T c 
this dependence is a universal function of u — (3kM(f3). This prediction is 
consistent with figure 2. 

The effect of pinning of the most important spins (Fig. 3) can be also 
described analytically and it follows that there occurs a discontinuous phase 
transition from the "spin up" phase to the "spin down" phase by a well 
defined critical number of pinned spins. 

Pinning one spin of degree k to the state S = — 1 is equivalent in MFA to 
lowering the mean magnetization by fj,{k)/N and introducing to the system 
the external magnetic field of the magnitude J = 1 that is oriented "down". 
This field is felt only by k other spins thus its mean value for the whole 
system equals b(k) = —k/N. If we pin j of the most connected spins then it 
means that we pin all spins of the degree k > k where 



It follows that we decrease the mean magnetization (per one spin) by 



M 




j = NA 




POO 

SM(j)=A (n(k) + l)k 
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where we have assumed that all pinned spins were completely ordered before 
pinning, fi(k) = 1. The effective internal field B(j) acting in the system from 
the pinned spins can be calculated as 

B(j) = -A j°° fi(k)k~ 2 dk ~ -2mJJjN . 

J K 



It is important to stress that 5M oc l/N but B(j) oc yl/N and this internal 
field causes much larger decrease of magnetization AM than the direct effect 
of pinning. For small j we have AM(j) = xB{j) where x is the system 
initial susceptibility. In MFA it can be calculated as the mean value of the 
derivative of local magnetization ji{k) over the field B(j) and it equals to 



2m 2 r 00 dk 

X = 



T Jm k 3 cosh 2 (f3M k) 

The last integral cannot be calculated analytically. Pinning several spins 
means influencing the system by a large internal field B(j) that can even 
cause a flip of the total magnetization. The value for the minimal flipping 
field B* can be obtained from the equation 

M* = A AT 3 tanh[/?(M*A; + B*)\dk 

J m 

and the corresponding tangency relation at the critical point M*: 

roo 

l=PA k- 2 cosh- 2 [p(M*k + B*)]dk . 

Jm 

The above equations form imply conditions for B* and M*. After some 
algebra one can find the following relation between B* and M*\ 

B* = -mM* +Tartanh(M*/2). 

Combining the last result with the relation for B(j) we get the following 
equation for the minimal number j* of pinned spins needed to invert the 
system magnetization: 

T M* 

— artanh— ~ (M*/2)(l - T/T c ) 

where the last approximation is valid only for the small M* = 0.03 ... 0.1 of 
Fig. 3 and we used the above T c = 2m. This prediction j*/N ~ 10~ 4 . . . 10~ 3 
agrees reasonably with our simulations. 
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In summary, we combined the Ising model with the Barabasi-Albert net- 
work and found that depending on the convincing power one either has a 
majority opinion or two equally widespread opposing opinions, i.e. a Curie 
point. 

We thank DAAD and the special grant of Warsaw University of Tech- 
nology to support the visit to Cologne University, the Julich supercomputer 
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Barabasi-Albert Ising model on 2 million nodes, m = 1 ,2,3,4,5 from left fo right; time = 1 00 or 500 



0.1 



0.01 r 



0.001 



0.0001 



° * x x + 
n * X 



□ X X 



1e+06 



5 10 15 20 25 30 35 40 

temperature 

Effective temperature, m=5, t=500, as a function of the number of nodes in the network 
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Figure 1: (a) M versus temperature for 2 million nodes and various m, for 
250 <t< 500 (shorter times t were used far below T c ). (For m = 1 even 60 
million nodes were simulated.) (b) Effective T c versus m + iV for m = 5 and 
various N, averaged over up to 1000 samples. 
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One network with N = 4000 at T = 2, 9, and 16 (M = 0.995, 0.509 and 0.147), and tanh(kM/T) 




kM/T 



Figure 2: Correlation between the number of neighbours and the local 
magnetization < Si > for one network of N — 4000 at T = 2, 9 and 16. 
Average over 250 < t < 500 iterations. The curve is the mean field prediction 
tanh((3kM). N = 2000000 gives similar effects. 
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Figure 3: Total magnetization versus time, summed over 100 networks of 
N = 30, 000 when after every 50 iterations the most-connected free spin is 
forced down permanently. For lower temperatures the sign change of the 
magnetization happens later. N = 2000000 gives similar effects. 
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